Finite element formulation for shear modulus reconstruction in transient elastography

In order to image the shear modulus in soft tissue, for medical diagnosis, given one component of measured displacements as a function of time on an imaging plane, two related direct finite element-based inversion algorithms are presented. One algorithm is based on the governing equations expressed in the frequency domain, and the other is in the time domain. The algorithms consider the complete equations of isotropic, small deformation, elasto-dynamics, where the hydrostatic stress is also treated as an unknown. The algorithms reconstruct both the shear modulus and hydrostatic stress fields, and regularization is used to stabilize the hydrostatic stress recovery. An algorithm is also developed for reconstructing the second displacement component, while simultaneous finding a smooth approximation to the measured displacement component to reduce noise. Shear modulus reconstruction results from both algorithms, using experimental ultrasound measurements on a tissue-mimicking phantom, are presented, and the merits and drawbacks of each algorithm are discussed.